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The positions of the Old Persian, Elamite, and Babylonian versions of the major trilingual inscription DB on the 
rock at Bisotün. Source: King and Thompson, pl. VI; corrected by Borger, fig. 2; adapted by R. Schmitt 
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2011. The Fukushima Daiichi Nuclear Disaster 





2011. The Fukushima Daiichi Nuclear Disaster 








2010, Deepwater Horizon Drilling Rig at Macondo Prospect 


Inextinguishable Fire, Sinking two days later, leaving the well gushing at the seabed. 








1986. Chernobyl Nuclear Disaster 


Ukrainian Soviet Socialist Republic, Soviet Union 





Reactor Safety Flaws and Human Actions 


Similar to Chernobyl Nuclear Power Plant 





2004, 2009 Closed two-unit RBMK-1500 Reactors 





Field Station Detail: 


Ignalina Nuclear Station Earthquake Early Warning 








The Yucca Mountain Nuclear Waste Repository 





The Yucca Mountain Nuclear Waste Repository 
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Hanford Radioactive Containment Site 





Hanford site represents two-thirds of the nation's high level radioactive waste by volume 
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Hanford Radioactive Containment Site 








Waste Isolation Pilot Plant, New Mexico, United States 





Radioactive Containment Casks arriving at WIPP 








Waste Isolation Pilot Plant 





Radioactive Waste is placed in Rooms Underground 
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Waste Isolation Pilot Plant, New Mexico, United States 
2014. Radioactive Leaks from Damaged Storage Drum 


Waste Explosion, Airborne Release of Radiological Material 








The Japanese government publishes а 
national seismic hazard map like this 
every year. But since 1979, earthquakes 











that have caused 10 or more fatalities in "Ри 
Japan have occurred in places it — -n 
designates low risk. ғаны 


1993 


78 (230) 





Fault plane 


011 Tohoku 


earthquake 
Magnitude-9.1 
(>27,000 dead 

or missing) 


1995 


7.3 (6,437) 
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Tonankai 
Tokai 
















San Andreas Fault System 


Slip segments during the 1906, 
1872 and 1857 earthquakes. 


Segments on which slip occurred 
during smaller earthquakes. 


Segments on which fault creep occurs. 








Italian Seismicity 
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3.2. The Cross-Covariance Function 


Time-distance helioseismology is based on the measurement of the cross-covariance between the 
Doppler signals at two points r, and r; on the solar surface, 


a 
Се, г.) = | dr ¢(r, Же, + D, (5) 
0 


where t is the correlation time lag. Figure ба shows а cross-covariance function measured from 
144 days of МОТ medium-degree data. The cross-covariance has been averaged over many pairs 
of points (r2, rı) and is presented as a function of the heliocentric angle between these two points. 
This diagram is known as the time-distance diagram. The cross-covariance is essentially a phase 
coherent average of the random oscillations (Bogdan 1997, and Supplemental Video 2). It is a 
solar seismogram: it provides a way to measure wave travel times between two surface locations. 
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Figure 4 Distribution of the епа nodes adopted for the tomogra] 
inversion in plan view (a) and vertical cross section (b). The open tr 
gles in (a) denote the four Apollo seismic stations. 











Figure9 Vertical cross sections of P-wave tomography along the 
four profiles shown on the map. Red and blue colors denote low and 
high velocities, respectively. The scale of velocity perturbations 
relative to the 1-D velocity model (Figure 3) is shown below (c). 
White dots denote moonquakes occurring within 150 km width of ~ \ ] / мл 
each profile. Open triangles on the map denote the four Apollo seis- “Т 7 60°8 
mic stations. — 
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The current velocity distribution is updated according 


I3s 





to the vector computed, and the process is repeated 
38 





until the Subsurface Map stabilizes. 


The ray path from a source to the station is defined by 
following the steepest gradient in the time grid from the 
source. From the ray geometry of all source-receiver pairs the -40' 
matrix describing the forward tomographic problem can be 
computed: 





г = GAu, (1) 


where г is the travel time residual vector, ie, the 45! 
difference between the observed travel times and the 
travel times computed in the current velocity distribution. 
Au is the slowness (1/velocity) perturbation vector. Each 
Gp; coefficient is the length of the ray p sampling cell j. — = A 
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Initial 3-2 Velocity Distribution: 
Starting Model designed using a 
subsurface simulation grid 

(nx, ny, nz) 


То image a region we use an estimation 
grid of cubic cells describing the scale 
we want to resolve (mx, my, mz) 


The current velocity distribution is updated according to the 
vector computed, and the process is repeated until the velocity 
model stabilizes. 





Depth [km] 


We give а quick overview of the theory of systems of linear, algebraic equations. This 
is not intended to serve as a text on linear algebra, merely a review of some important 
concepts. Detailed discussions of various aspect of this material can be found in [43] or 
[78]. 


Consider a system of m equations for n unknowns, (£1,..., £n): 


а,171 + а1212+ +++ cb Gigs =Y 
а2171 + а22т2+ Famn = yo 

| | (1.16) 
Qm1 11 + (45232 + sco Ст = Ym- 


There аге four questions which require answers: 


4.2.1 Inverse problem 


The inverse problem corresponding to the forward problem given 
by 2 is nonlinear. We define the a posteriori probability density 


function (PDF) as 
ом(т)= 
ог 2 
1 1—1 (se <) N—1 5, (% =) 
(EESE) ,уа(5—: 
N-1 pio 2 — M-1 пог 2 
S; (lı - 1, Ak — А 
+ ( +> ( Е ) (7) 
j=0 k=0 


Existence: 
For a given m-vector (уј,..., Ут) does there exist an n-vector (71,...,2,) which 
satisfies the equations in (1.16)? 


Uniqueness: 
When а solution exists is it unique? More generally, describe the space of solutions. 


Solve in practice: 
Give an algorithm to find approximations for the solutions of (1.16) and criteria 
to select a solution when there is more than one. 


Stability: 
How sensitive is the solution to small variations in the coefficients (a;;) or the right 
hand side (y;)? 


It is а somewhat unexpected, but very important fact that these issues are in practice, 
rather independent of one another. 
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